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The MINOS Collaboration continues its search for u e appearance in the NuMI (Neutrinos at the Main Injector) 
beam at Fermilab. Neutrinos in the beam interact in the Near Detector, located 1 km from the beam source, 
allowing us to characterize the backgrounds present in our analysis. In particular, we can estimate the number 
of v e candidate events we expect to see in the Far Detector (735 km away, in the Soudan mine in northern 
Minnesota) in the presence or absence of f M — > u e oscillation. Recent efforts to improve the sensitivity of the 
analysis, including upgrades to the event identification algorithm and fitting procedure, are discussed, and the 
latest results from the search are presented. 



1. Introduction and Motivation 

Neutrino oscillation is now an established phenomenon, thanks to a variety of experiments [IHZ] which have 
probed various apparent anomalies in expected neutrino spectra from the sun, atmosphere, nuclear reactors, 
and accelerator-based beams. The three weak eigenstates - u e , v^, v T — are related to the three mass eigenstates 
- v\, i>2, vj, - via the PMNS (for Pontecorvo, Maki, Nakagawa, and Sakata) matrix, a unitary mixing matrix U 
parametrized by three mixing angles - f?i 2 , #13, 6*23 - and three CP-violating phases -8, a±, a 2 [8j. No hints 
exist yet as to the value of 8, and oscillation experiments are not sensitive to the latter two CP-violating phases, 
which are nonzero only if neutrinos are their own antiparticles (i.e., Majorana particles). So far, two of the 
mixing angles, 0i 2 and 23 , have been measured [H 01 E] and are large. The value of #13, on the other hand, 
appears to be quite small. It was first significantly constrained by the CHOOZ experiment in 1999 [9|; the 
result is shown in Fig. [T] These proceedings report the latest improvements on this constraint by the MINOS 
experiment. 

The probability of a neutrino of energy E created with flavor a being detected with flavor /3 after traveling 
some distance L in vacuum - to wit, the oscillation probability - is given by: 



P{v a -> vp) = 



where rrij is the mass of the j-th neutrino mass eigenstate. The final expression is found to depend not on the 
absolute squared masses but on two independent squared-mass differences - Am^j and Am§ 2 . While both the 
sign and magnitude of the first of these are known [I] , only the magnitude of the second squared-mass difference 
is known [B]. That is, we do not know if 1/3 is more or less massive than v\ and v-i\ the first possibility is referred 
to as the normal mass hierarchy and the latter, the inverted mass hierarchy. 
The oscillation mode of present interest is v e . To leading order: 



P{y^ -> u e ) « sin 2 (26»i 3 ) sin 2 (6> 23 ) sin 2 (1 .27 Aml 2 L / E) , (2) 

with Am§ 2 in units of eV 2 , L in km, and E in GeV. When we account for the propagation of the neutrinos 
through matter, additional, higher order terms appear in Eq. [2] [TU]. The full expression depends on the CP- 
violating phase 8 and the mass hierarchy; accordingly, the results presented here will be sensitive to these 
parameters. 

Under certain conditions leptonic CP violation can lead to non-conservation of baryon number, which 
in turn can help explain the observed matter-antimatter asymmetry in our universe. Access to the value of the 
leptonic CP-violating phase, namely 8, is controlled by the size of 6*13. Should #13 be sufficiently large, one may 
be able to measure 8 with current or next-generation long-baseline neutrino oscillation experiments. Finding, 
on the other hand, that #13 is zero could point to a new symmetry. 
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Figure 1: Exclusion contours from the CHOOZ experiment constraining the value of the #13 mixing angle (horizontal 
axis) as a function of Am| 2 (vertical axis). Taken from Ref. [5]. 



2. The MINOS Experiment 



The MINOS experiment consists of three key parts: the Neutrinos at the Main Injector (NuMI) neutrino 
beam at Fermilab, the Near Detector (ND) at Fermilab, and the Far Detector (FD) in Soudan, MN. The 
NuMI beam [12] is produced by striking a graphite target with 120-GeV protons from Fermilab's Main Injector, 
focusing the resulting secondary hadrons (mostly pions, some kaons) using pulsed-current magnetic horns, and 
allowing these hadrons to decay leptonically in flight. Although the horns' polarities are set to focus 7T + and 
K + and defocus tt~ and K~ , a small fraction of negative hadrons travel along the horns' axis, avoid being 
defocused, and remain in the beam. The primary decay mode for these hadrons is or [i~v^\ a very small 

fraction decay to e + v e or e~T> e . Ultimately, the resulting neutrino beam is 92% and 7% v^; v e and v e make 
up the remainder. 

The MINOS detectors are functionally identical, magnetized tracking calorimeters. The 980-ton ND is located 
approximately 1 km downstream of the NuMI target and measures the neutrino beam's initial composition and 
energy spectrum. The 5.4-kton FD is located 735 km downstream in the Soudan mine in northern Minnesota 
and measures the oscillated beam's composition and energy spectrum. Both detectors consist of alternating 
layers of 2.54-cm thick steel and 1-cm thick plastic scintillator. The scintillator is segmented into optically- 
isolated, 4.1-cm wide strips and serves as the active portion of each detector. Light from the strips is read out 
via optical fibers and multi-anode photomultiplier tubes. Further details can be found in Ref. |13j . 
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3. Neutrino Interactions in the MINOS Detectors 
3.1 . Event Types 

Though a multi-purpose experiment, MINOS was designed and optimized to study disappearance. The 
event type of interest in such an analysis is a charged-current (CC) interaction: 

+ N -> p- + X, (3) 

where N represents a nucleon or a quark inside a nucleon and X represents the resulting hadronic recoil system. 
The track produced by the muon typically extends far beyond the shower activity of the recoil system (see left 
of Fig. |2j, enabling straightforward identification of such an event. In the analysis at hand, however, the event 
type of interest is a v e CC interaction: 

u e + N -> e~ + X. (4) 

In the MINOS detectors, the electromagnetic shower generated by the electron has a radiation length of 1.76 cm 
- ~30% less than the thickness of a steel layer - and a Moliere radius of 3.7 cm - ~10% less than the width of a 
scintillator strip. As a result, the electromagnetic shower becomes fully interspersed with the hadronic shower 
(see center of Fig. [2| . The challenge, then, is to distinguish this type of event from the primary background 
event type - the neutral-current (NC) interaction (see right of Fig. [2]): 

v x + N -> u x + X. (5) 



Event II », Event | I NC 
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Figure 2: Displays of simulated events in the MINOS detectors. Pulse heights shown are calibrated such that 1 unit is 
the amount of energy deposited by a minimum ionizing particle. 



3.2. Event Identification 

This is the third v e appearance search conducted by MINOS. In all three analyses, event selection consisted 
of three general stages. First, we required that the event took place within a fiducial volume. Next, we rejected 
events that were clearly or likely background. This "pre-selection" eliminated most CC events by requiring 
that any reconstructed track in the event be no more than 24 planes long and extend no more than 15 planes 
beyond the hadronic shower. We also required that an event have at least five contiguous planes with energy 
greater than half that deposited on average by a minimum ionizing particle. Events were required to be in-time 
with and to point in a direction consistent with the neutrino beam pulses. The reconstructed energy of the 
event had to lie between 1 and 8 GeV, as most events below 1 GeV are from NC interactions and P{v^ — > v e ) 
is negligibly small above 8 GeV. 

The last stage of the selection involved the use of a sophisticated event identification technique to separate 
signal from background. In our first two v e appearance analyses 115) . we used a neural network (named 
Artificial Neural Network, or ANN) whose input variables were reconstructed quantities characterizing the 
longitudinal and transverse energy deposition profiles of an event. The present analysis utilizes instead a 
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novel technique named Library Event Matching (LEM) that uses raw energy deposition information instead 
of reconstructed quantities and is based on a pattern recognition algorithm. We begin by comparing a given 
candidate event to each of 5 x 10 7 simulated events (the "library"), of which ~ 3 x 10 7 are background events 
and ~2x 10 7 are signal events. Intuitively, the comparison determines the extent to which the topologies and 
strip pulse heights of the candidate and library events are similar. Quantitatively, this is performed using a 
likelihood: 



Nstrips N planc3 

log£= lo s 

3=1 fc=l 

where P(n, A) is the Poisson probability for observing n given mean A and is the charge (in photoelectrons) 
deposited in strip j of plane k of event x (candidate or library). Once this matching process is completed for 
a given candidate event, the library events are ranked by their log£ values. The top 50 of these, designated 
the "best matches," are summarized by three variables: (1) the fraction that are signal (v e CC) events, (2) for 
these signal events, the average fraction of total event energy in the hadronic recoil system, and (3) for these 
signal events, the average fraction of deposited charge that overlaps between each of them and the candidate 
event. These variables, along with the reconstructed energy of the candidate event, are given as inputs to a 
neural network; the output is the LEM discriminant, shown in Fig. [3] 



(6) 




Figure 3: LEM discriminant distributions for signal and background events passing the fiducial and pre-selection cuts. 



4. The Analysis 
4.1 . Overview 

The present analysis represents three improvements over the most recent previous analysis, the effects of which 
are summarized by Fig. |4j An additional exposure of 1.2 x 10 20 protons on the NuMI target provides a modest 
17% increase in the size of our data set. Using LEM instead of ANN provides a gain in sensitivity to 6*13 = of 
about 15%. Finally, instead of cutting at the optimal value of the discriminant (found to be 0.7 for LEM) and 
comparing the number of events above that value to the expectation (as we did in the previous two analyses) , we 
perform a fit to the LEM discriminant distribution (optimal binning was found to be 0.6-0.7, 0.7-0.8, 0.8-1.0), 
obtaining an additional 12% increase in sensitivity. The basic approach to making a background prediction, 
however, remains unchanged - count the events observed in the ND and use the FD and ND simulations to 
"extrapolate" this to the FD. 



4.2. Background Prediction 



Events due to NC interactions constitute the primary background in this analysis. Contributions to the 
background from CC and intrinsic beam v e CC interactions are smaller but non-negligible; the — > v T 
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Figure 4: Projected 90% C.L. sensitivities to #13 = assuming the normal mass hierarchy. The contour labeled "ANN 
> 0.7" is the sensitivity achieved by repeating the last analysis with the additional data included. The "LEM > 0.7" 
contour is the sensitivity for a counting experiment performed with LEM instead of ANN using the same full data set. 
Finally, the "LEM 3 bins" sensitivity is obtained by performing a shape fit to the LEM distribution. 




oscillation mode allows the possibility of v T CC interactions, an additional source of background. Such a multi- 
component background complicates the aforementioned extrapolation process, as each of the components is 
affected differently by the transition from ND to FD (e.g., NC interaction rates are not affected by oscillations 
while the various CC interaction rates are affected to differing extents). Consequently, we perform the back- 
ground prediction componcnt-by-componcnt. In particular, the background prediction for component a (NC, 
Ufj, CC, v e CC) in analysis bin i is given by: 



171 ^Simulation 

FD Pred t cted = jy pData «,i 

a > 1 a < 1 at ijSirnulation ' 



a. i 



where the individual FD and ND terms stand for event counts. In a manner similar to that of the v r CC 
interactions, signal v e CC interactions "appear" due to oscillations (in this case, via the — ¥ v e mode). As 
such, these event types appear only in the FD and cannot be extrapolated from the ND using Eq. [7J we use 
a more involved but conceptually similar approach based on the observed ND CC spectrum to predict the 
event counts due to these interactions. 

The ND^ ta term stands out on the right-hand-side of Eq. [7] While the total number of events in the ND 
data is straightforwardly obtained (simply apply the relevant cuts and count them) , the means of determining 
the individual contributions of the three background event types to this number are far from obvious. Our 
solution to this problem takes advantage of the flexibility of the NuMI beam configuration. In addition to 
taking data in the standard configuration, we also take data (albeit with considerably smaller exposures) in 
the so-called "horn-off" (focusing horns turned off) and "high-energy" (target pulled farther away from the 
focusing horns, horn currents increased slightly) configurations. These configurations produce neutrino beams 
with different true energy spectra (see Fig. [5]), which in turn produce activity in the ND with different relative 
fractions of NC, CC, and intrinsic beam v e CC events. 

Combining the total ND data rate in each of the three configurations with the simulated background compo- 
nent ratios between configurations, we perform a fit and obtain the ND®f ta , the data background components 
in the standard configuration, as shown in Fig. [6j 

The FD-to-ND ratio in Eq. [7] is taken directly from the simulation and is responsible for the extrapolation 
from ND to FD. It corrects for a variety of effects, including flux differences (1/i? 2 fall-off, acceptance, decay 
kinematics, focusing), fiducial volume effects, energy smearing, oscillations, and detector effects. Similarities 
between the detectors, including associated similarities between their simulations, allow much of the systematic 
uncertainty to cancel to first order in this ratio. The resulting values of this ratio, with systematic error bars, 
are shown in Fig. [7] for each of the three ND-based background component types, as a function of reconstructed 
energy. 

With this information, we are now able to make the background and signal predictions. In particular, we 
have the ability to adjust the predictions for different sets of oscillation parameters. To illustrate, Table [I] lists 
the predicted counts of events with LEM > 0.7 for two reference values of #13, assuming a commonly used set 
of values for other oscillation parameters. 
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Figure 5: Exposure-normalized energy spectra of all events occurring in the ND fiducial volume, for each of three different 
beam configurations described above. 
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Figure 6: The energy spectrum of all events in the ND used in the official analysis (i.e., those with LEM > 0.6 and 
passing the fiducial and pre-selection cuts). The NC, CC, and CC spectra are obtained by the technique described 
above. 



4.3. The FD Data 



After applying the fiducial and pre-selection cuts to the events in the FD data, we obtain the energy spectrum 
shown in Fig. [8] Although the pre-selection cuts events with reconstructed energy above 8 GeV, we performed 
the prediction out to 12 GeV and compared it to the data to test the robustness of our background prediction 
method. 

In the LEM discriminant distribution, we expect negligible signal in the LEM < 0.5 region and thus can 
use it as a test of the entire analysis chain. We predict 370 ± 19(stat.) background events with LEM < 0.5 
and observe 377 - well within statistical error. Sufficiently assured of the soundness of our analysis method, 
we proceed to unblind the signal-enhanced region of the LEM discriminant distribution. The full spectrum is 
shown in Fig. [9] For reference only, we can count the number of events with a LEM discriminant above the 
optimal cut of 0.7; we predict 49.5 ± 2.8(syst.) ± 7.0(stat.) background events and observe 62 events, an excess 
of approximately 1.1a. 
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Figure 7: The FD-to-ND ratio for the three background event classes as a function of reconstructed energy. These ratios 
were calculated for events with LEM > 0.7 (the optimal single cut on the discriminant) and are thus for illustration 
only; the official analysis uses ratios calculated in each bin of the LEM discriminant. The "Run Periods" referred to in 
the plot are periods of experiment running with slightly different conditions; these conditions are accounted for in the 
simulation, thus explaining why the values of the ratio differ slightly between Run Periods. 



Table I: Predicted event counts for events with LEM > 0.7. These numbers were calculated assuming 8 — 0, 6*23 = tt/4, 
and Am| 2 = 2.32 x 10~ 3 eV 2 (i.e., normal mass hierarchy). The line labeled "v e CC" refers to the intrinsic beam v e 
component. The nonzero signal in the #13 = case is due to AmJi-driven oscillations. 
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4.4. Fitting to the Data 

The official fit to the data is performed in three LEM bins and five reconstructed energy bins. As an example, 
Fig. [TO] shows the observed reconstructed energy spectra in the three LEM bins superimposed on the best-fit 
signal and predicted background under the conditions indicated. Our official result, described in Ref. |16| . 



is given in Fig. 11 These contours were calculated according to the Feldman-Cousins prescription [17] and 
thus provide correct coverage near the physical boundary at # 13 = 0. In addition, uncertainties in the other 
oscillation parameters have been accounted for in generating these contours. Assuming 5 = and 6*23 = tt/4, we 
find that sin 2 (26*13) = 0.04llg'g 3 i (Q.Q79j^"gg3) for the normal (inverted) mass hierarchy. Under the same set 
of assumptions, sin 2 (26*13) < 0.12(0.20) at 90% C.L. for the normal (inverted) mass hierarchy. Finally, these 
results exclude the null hypothesis of sin 2 (26*13) = at 89% C.L. 
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Figure 8: Predicted (assuming #13 = 0) and observed reconstructed energy spectra in the FD. 
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Figure 9: Predicted (assuming #13 = 0) and observed LEM discriminant distributions in the FD. 



References 

1 B. Aharmim et al. (SNO), Phys. Rev. Lett. 101, 111301 (2008). 

2 T. Araki et al. (KamLAND), Phys. Rev. Lett. 94, 081801 (2005). 

3 Y. Ashie et al. (Super-Kamiokande), Phys. Rev. Lett. 93, 101801 (2004). 

4 Y. Ashie et al. (Super-Kamiokande), Phys. Rev. D 71, 112005 (2005). 

5 K. Abe et al. (Super-Kamiokande), Phys. Rev. Lett. 97, 171801 (2006). 

6 P. Adamson et al. (MINOS), Phys. Rev. Lett. 106, 181801 (2011). 

7 P. Adamson et al. (MINOS), Phys. Rev. Lett. 107, 021801 (2011). 

8 B. Pontecorvo, JETP 34, 172 (1958); V. N. Gribov and B. Pontecorvo, Phys. Lett. B 28, 493 (1969); Z. 
Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

9 M. Apollonio et al. (CHOOZ), Phys. Lett. B 466, 415 (1999). 

10 E. K. Akhmedov et al., J. High En. Phys. 04, 078 (2004). 

11 A. D. Sakharov, JETP Lett. 5, 24 (1967). 

12 K. Anderson et al., FERMILAB-DESIGN-1998-01 (1998). 

13 D.G. Michael et al., Nucl. Inst. & Meth. A 596, 190 (2008). 

14 P. Adamson et al. (MINOS), Phys. Rev. Lett. 103, 261802 (2009). 

15 P. Adamson et al. (MINOS), Phys. Rev. D 82, 051102 (2010). 



16 P. Adamson et al. (MINOS), arXiv:1108.0015 vl [hep-ex] (2011). Accepted by Phys. Rev. Lett 



17 G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998). 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



9 



MINOS Far Detector Best Fit (0.6 < LEM < 0.7) 
-m- FD Data 

Background 

Signal 

sin a (26 13 )=0.041 , Am^>0, S cp =0 
Merged for Fit 




2 3 4 5 6 7 
Reconstructed Energy (GeV) 



t— 20 
O 

Q. 

ii 5 ; 

X 
CM 

S? 10 : 

o • 

c - 

0) 



MINOS Far Detector Best Fit (0.7 < LEM < 0.8) 
-■- FD Data 

Background 

; Signal 

Hi IH sin a (2e.,)=0.041 , d<>0, Sr. p =0 



Merged for Fit 



2 3 4 5 6 7 
Reconstructed Energy (GeV) 



O 

Q. 

*>1 



;io 



MINOS Far Detector Best Fit (LEM > 0.8) 
-■- FD Data 

Background 

Signal 

sin 2 (26 13 )=0.041, Arr^ 2 >0, S cp =0 Z 
Merged for Fit 




2 3 4 5 6 7 
Reconstructed Energy (GeV) 



Figure 10: Stacked background prediction and best-fit signal assuming 8 = 0, 023 = 7t/4, and Ami 2 = 2.32 x 10" 3 eV 2 
with data superimposed. Note that each reconstructed energy spectrum is for a different bin of the LEM discriminant. 




Figure 11: Allowed ranges and best fits for 2 sin 2 (#23) sin 2 (2f5i3) as a function of S. The upper (lower) panel assumes 
the normal (inverted) mass hierarchy. The CHOOZ 90% C.L. upper limit shown assumes 623 = 7r/4 and Am^ = 
2.32 x 10" 3 eV 2 . 



